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<3 Abstract 
O 

Observations of high redshift supernovae imply an accelerating Universe 
which can only be explained by an unusual energy component such as vacuum 
energy or quintessence. To assess the ability of current and future supernova 
data to constrain the properties of the dark energy, we allow its density to 
have arbitrary time-dependence, px{ z )- This leads to an equation of state for 
the dark energy, Wx(z) = px(z)/ px(z), which is a free function of redshift 
\ z. We find that current type la supernova (SNe la) data are consistent with 

a cosmological constant, with large uncertainties at z > 0.5. We show that 
Px{z)l Px(z = 0) can be measured reasonably well to about z = 1.5 using type 
la supernova data from realistic future SN la pencil beam surveys, provided 
that the weak energy condition (energy density of matter is nonnegative for any 
observer) is imposed. 

While it is only possible to differentiate between different models (say, 
quintessence and k-essence) at z < 1.5 using realistic data, the correct trend in 
the time-dependence of the dark energy density can be clearly detected out to 
z — 2, even in the presence of plausible systematic effects. This would allow us 
to determine whether the dark energy is a cosmological constant, or some exotic 
form of energy with a time- dependent density. 
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1. Introduction 



Most of the present energy content of our universe is unknown QBahcall et al. 1999|) . 
Distance-redshift relations derived from cosmological standard candles at redshifts between 
zero and a few are the most sensitive probes of the equation of state of the universe, which 
allows us to constrain the energy content of the universe. 

Type la supernovae (SNe la) are our best candidates for cosmological standard candles. 
They can be calibrated to have small dispersion in their intrinsic luminosities (|Phillips| 
|1993| , [Riess, Press, fc Kirshner 1995| ). The data from two independent observational teams, 
the High-z SN Search (Schmidt et al.) and the Supernova Cosmology Project (Perlmutter et 
al.), seem to suggest that our universe has a significant vacuum energy content ( |Garnavich] 
|et al. 1998| , [Riess et al. T998j Perlmutter et al. 1999Q . 



While a cosmological constant term (vacuum energy) in the Einstein equations 
provides the simplest explanation of the current SN la data, other forms of energy 
(quintessence, dark energy, etc) have also been studied (e.g. [White 19981 , IGarnavich et al 



[TOj Bernhardt, Wang, fc Zlatev 1999|, |Kfstathiou 1999; [Podariu fc Ratra 2000[ |Waga fc 



Frieman 2000| ) and are consistent with current data as well. Since cosmology has matured 
into a phenomenological science at the turn of the new millennium, observational data will 
dominate aesthetics in the selection of cosmological models. 

So far, most cosmologists have assumed time-independent equations of state, i.e., 
power-law dark energy density (see §2), in the context of constraining the energy content 
of the universe. In this paper, we allow the dark energy density of the universe to be a 
free function of redshift, i.e., arbitrary equation of state. We derive constraints on the 
dark energy density from current SN data, and assess future prospects of measuring the 
dark energy density using simulated data from realistic future SN la surveys. We expect a 
model-independent measurement of the dark energy density as a function of redshift to be 
very useful in our quest of the unknown energy contents of the universe. 



2. Parametrization of the dark energy density 



In a smooth Friedmann-Robertson- Walker (FRW) universe, the metric is given by 
ds 2 = dt 2 — a 2 (t)[dr 2 /(l — kr 2 ) + r 2 (d6 2 + sin 2 9d(j) 2 )}, where a(t) is the cosmic scale factor, 
and k is the global curvature parameter. The cosmological redshift z is given by 1 + z = 1/a. 

To make model-independent measurements of the equation of state, we replace the 
vacuum energy density p\ with p x = p x f(z) in the total matter density of the universe: 



P(z) = Pl{l + zf + pl{l + z) 2 + p x f{z) 

^ m (l + zY + Sl k (l + z) 2 + Sl x f(z) 



(i; 
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where the superscript "0" indicates present values, f(z = 0) = 1, and 
Q k = 1 — Q m — Qx — ~k/HQ. If the unknown energy is due to a cosmological 
constant A, f(z) = 1. Clearly, the function f(z) is a very good probe of the nature of the 
unknown energy. 



The comoving distance r is given by flWeinberg 19721 ) 



r(z) = cH, 



K ' 



r(z; n m , n x , f) = / dz—— 

Jo tejKZ') 



E(z) 



Q m {i + z'Y + n x f{ z ') + n k {i + z 



l\2 



1/2 



(2) 

(3) 
(4) 



where 



«S'(x) = sinh(x) 

= ^, 

= sin fx), 



Q k > 
< 0. 



(5) 



The angular diameter distance is given by dj^iz) = r(z)/(l + z), and the luminosity distance 
is given by d L (z) = (1 + z) 2 d A {z). 



Einstein's equations in a FRW metric, together with the first law of thermodynamics 
give us 

rln 

(6) 



dp 



For unknown energy px{z), we find 

Px{z) = p° x f(z) 
Px(z) 

Now we can write the equation of state of the unknown energy as 



l -{l + z)f\z)-f{z) 



(7) 



/ \ _ Px(z) 1 f (z) 

w{z) = — — = -(1 + z)—— - 1. (8 
Px\z) 3 f(z) 

A constant equation of state corresponds to f(z) oc (1 + z) a , where a is a constant. The 
values a = 0, a = 3, and a = 4 correspond to a cosmological constant, matter, and 
radiation respectively. 

To obtain accelerated expansion, we need p + 3p < 0. Since p x + Spx = 
p° x [(1 + z)f - 2/], this implies a< 2 for f(z) oc (1 + z) a . 
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The weak energy condition states that for all physically reasonable classical matter, 
the energy density of matter as measured by any observer is nonnegative ( |Wald 1984|) . For 
a perfect fluid, the weak energy condition will be satisfied if and only if 



This leads to (see Eq. [ 



p + p > 0. 



/'(*) > o. 



(9) 



(10) 



The weak energy condition imposes strong constraints on the jointly estimated cosmological 
parameters Q m , Oa, and the dark energy density f(z). 

Note that the co moving distance r(z) depends on the equation of state of X through 



T(z) 
Hence 



dz' 



n m {i + z'y + n x f(z') + n k (i + z 



l\2 



-1/2 



H 



dr 



o y/l — kr z 



(11) 



v'{z)= n m (i + z) 3 + n x f(z) + n k (i + z) 



-1/2 



Hn dr , , 

mi:- ( 12 ) 



[l + fi fe # V] 1/2 dz 



To measure f(z) directly from data, we need to evaluate the derivative of the distance r(z) 
with respect to redshift z. 

To avoid taking derivatives of noisy data, we can parametrize f(z) with its values at n 
equally spaced redshifts Zi, and assume that f(z) is given by linear interpolations at other 
values of z. We write 



Zi - z 



\ Zi Z; L _\ J 

zq 0, z n z maX j 



h-i + 



Z — Zi-i 



fi ■ 



y Zi ^i— 1 / 
fo =1, fn = fn-1 



Z%-\ < Z < Zi 



(13) 



where fi (i = 1, 2, n — 1) are independent variables to be estimated from data. 



3. Parameter estimation 

The measured distance modulus for a SN la is 

$=H? + eW (14) 

where jJp is the theoretical prediction 



(15) 
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and is the uncertainty in the measurement, including observational errors and intrinsic 
scatters in the SN la absolute magnitudes. 

Denoting all the parameters to be fitted as s, we can estimate s using a x 2 statistic, 
with (Riess et al. 1998) 



2 



J l) ] 

x 2 (s) = e f "r •/•''• - e rp v "": 2 J , (is) 



°Mo,Z + a mz,l I °l 

where <t w is the estimated measurement error of the distance modulus, and a mz is the 
dispersion in the distance modulus due to the dispersion in galaxy redshift, a z , due to 
peculiar velocities and uncertainty in the galaxy redshift (for the Perlmutter et al. data, 
the dispersion due to peculiar velocities is included in cr m eff, i.e., <7 M0 ). Since 

5 (lddA 

^hTioU^T 2 ' (17) 

a rnz depends on the parameters s. The probability density function (PDF) for the 
parameters s is 

p(s) oc exp (-^ . (18) 

The normalized PDF is obtained by dividing the above expression by its sum over all 
possible values of the parameters s. 

In order to impose the weak energy condition f'(z) > 0, we compute the PDFs on 
a iV-dimensional grid for iV parameters. The PDF of a given parameter Sj is obtained 
by integrating over all possible values of the other N — 1 parameters. To reduce the 
computation time, we can integrate over the Hubble constant H analytically, and define a 
modified x 2 statistic, with 

~2 _ 2 Ci 



where 



X =X*-^(Ci + -lnlO), (19) 



where 



x* = i^ — yVp Vb) 

i a i 

i 

C2 = E4 (20) 
//* = ^ p (/i = h*) = 42.384 - 5 log h* + 5 log [H r(l + z)] . (21) 
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We take h* = 0.65. Our results are independent of the choice of h*. 

After experimenting with a number of different techniques, we developed an adaptive 
iteration method of estimating f(z) based on the requirement that f'(z) > (i.e., the 
weak energy condition is satisfied). Starting with the initial guess of f(z) — f(z — 0) = 1 
(a cosmological constant), we iteratively build up f(z) as parametrized by Eq.(|i~3"l) while 
minimizing y 2 . 



4. Constraints of px(z) from current SNe la data 



Wang (2000b) has combined the data of the High-z SN Search team QSchmidt et al.| 
|1998| , |Garnavich et al. 1998| , [Riess et al. 1998] ) and the Supernova Cosmology Project 
( Perlmutter et al. 1999| ), yielding a total of 92 SNe la. Using the method described in the 
previous section, we estimate Q m , Qx, and f(z) simultaneously by minimizing the modified 



X 



(see Eq.©). 



Fig.l shows the dimensionless dark energy density f(z) (as parametrized by Eq. (fT3|) ) 
measured from this set of 92 SNe la. It is consistent with f(z) = l (a cosmological constant), 
with large uncertainty beyond z > 0.5. Note that none of the error bars extend beneath 
f(z) = 1, because we have imposed the weak energy condition, i.e., f'(z) > 0, which implies 
that f(z) > f(z = 0) = 1. The estimated values of VL m = .3 (0, .9); tt x = 1.7 (.4, 2.2) are 



consistent with previous results QGarnavich et al. 1998| , [Riess et al. 1998] , Perlmutter et 



al. 1999| , |Wang 2000 b| ). The errors are estimated from the ranges of parameters for which 

y 2 = y 2 . +1. 

A A.mm ' 

Wang (2000b) found that when fit to a model with a cosmological constant as the dark 
energy, flux-averaging changes the best fit model to this data set of 92 SNe la. Without 
flux-averaging, the best fit model is a closed universe with fl m = 0.7 ± 0.4, and a vacuum 
energy density fraction of fl\ = 1.2 ± 0.5, consistent with the estimated values in Fig.l. 
The flux-averaged data yield fl m = 0.3 ± 0.6, and = 0.7 ± 0.7. This difference may have 
resulted from the data containing large redshift dependent uncertainties, which would have 
caused the results from the unbinned data to be biased. The effect of flux-averaging on the 
best fit model assuming an arbitrary dimensionless dark energy density f(z) will be studied 
elsewhere. 



Future cosmic microwave background (CMB) space missions MAP ( Pennett et al 



|1997[ ) and Planck ( pe Zotti et al. 1999| ), together with the galaxy redshift surveys SDSS 
( Cumi 1999|) and 2df (Palton et al. 2000]) , will give us exquisitely accurate measurements 
of the geometry of the universe and the matter density in the universe ([Eisenstein, Hu, fc| 
legmark 1999| , [lurner fc Tyson 1999| , |Wang, Spergel, fc Strauss 1999| ). SN data can provide 



the unique probe on the nature of dark energy by allowing us to measure how the dark 
energy density varies with time. 



- 7- 



Current cosmic microwave background (CMB) anisotropy measurements seem to 
indicate that we live in a flat universe (|de Bernardis et al. 2000| , |Balbi et al. 2000"!) . Cluster 
abundances strongly suggest a low matter density universe ( |Bahcall, Lubin, fc Dormari 
19951 ICarlberg et al. 1996| , Bahcall fc Fan 1998Q . fi m = 0.3 and fi A = 0.7 is the best fit 
model to current observational data. We will use Q m = 0.3 and Qx = 0.7 for our simulated 
data in the rest of this paper. 



5. Measuring px{z) from future SNe la data 



A large number of SNe la at z > 1 is critical in resolving the important systematic 
uncertainties of SNe la as cosmological standard candles, such as dust ( [Aguirrc 1999 ), 
gravitational lensing ( [Kantowski, Vaughan, fc Branch 1995| , [Wambsganss et al. 1997| , |Holz 
1998| IMetcalf fc Silk 1999| , |Wang 1999| , [Barber et al. 2000| ), and luminosity evolution ( [DreIT 



|Loredo, fc Wasserman 200~0| , [Riess et al. 1999| , [Wang 2000b|) , and in making SNe la useful 
probes of the dark energy content of the universe. The most efficient method of obtaining 
a large number of SNe la at z > 1 is conducting a supernova pencil beam survey on a 



dedicated large aperture telescope with a square degree field of view ( |Wang 2000ap . 



To study how well SN data can probe the dark energy density, let us consider two 
hypothetical dimensionless dark energy densities f q (z) and fk{z), given by 



/,(*) 
/*(*) 



,1.52 



(1 + z) 1 - 5 ' 
exp[0.9(l - e 



w g (z) = -1 + 0.5z 
)], w k (z) = 0.3(1 + z)e 



- 1 



(22) 



We have chosen f q (z) and fk(z) to represent quintessence (dw q /dz > 0) and k-essence 
(dwk/dz < 0) models respectively flCaldwell, Dave, fc Steinhardt 1998| , |Armendar iz- P icon 
[Mukhanov, fc Steinhardt 2000|) . Note that f q (z) and fk(z) satisfy the weak energy condition 
f'(z) > 0; they give an accelerating universe for z < 1.33 and z < 2 respectively. 

A feasible SN pencil beam survey (either from ground^ or from space), with a square 
degree field of view and for an effective observational period of one year, can yield almost 
2000 SNe la out to z = 2 ( |Wang 2000a| ). Let us combine the data from the SN pencil beam 
survey with SN data at smaller redshifts, so that there are a minimum of 50 SNe la per 0.1 
redshift interval at any redshift. This yields a total of 1966 SNe la for the quintessence and 
1898 SNe la for the k-essence model, up to z = 2 and for Q m = 0.3, Qx = 0.7. We simulate 
the data by placing perfect standard candles at random redshifts, with the total number 
per 0.1 redshift interval given as described above. Then we add intrinsic and observational 
dispersions which are Gaussian with zero mean and a variance of 0.20 magnitudes. A 



2 The SNe la at z > 1.5 will likely require follow up spectroscopy from space. 
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systematic shift in /i of dm sys z magnitudes is also added to mimic possible systematic 
errors as one goes to larger redshifts. 

Realistic SN la data should contain gravitational lensing noise. Wang (2000a,b) has 
shown that flux averaging is important in reducing the bias due to lensing. The effect of 
gravitational lensing and of flux averaging in the context of a general equation of state will 
be investigated elsewhere. 

Fig.2 illustrates how well we can recover the dark energy densities f q (z) and fk{z) in 
the absence of lensing noise and systematic shifts, when we apply our adaptive iteration 
method to 100 random data sets with a realistic dispersion of 0.2 magnitudes. To study 
the dependence of our results on the parametrization of f(z) [see Eq.(|l3"Dl, we show results 
for (a) n = 10, and (b) n = 6. The thick and thin solid lines are the assumed true f q (z) 
and fk(z) respectively. We have assumed that we know fl m + Qx — 1- The mean and 1-a 
errors of f(z) and Q m are estimated from averaging over the 100 random samples. The 
error of estimated Q m reflects the resolution of the program. For a given data set, the 
recovered f(z) and Q m are expected to fall within the errors with 68% probability. Clearly, 
the quintessence and k-essence models can be differentiated marginally for z < 1.5. For 
z > 1.5, the errors increase significantly, while the estimates become more biased, making it 
impossible to differentiate between the two models. However, the correct trend in the time 
variation of the dark energy density can be clearly detected out to z = 2. 

Fig.2 shows that the parametrization of f(z) with n = 10 yields less biased estimates 
than n = 6. This is as expected, since f(z) is more accurately parametrized as one increases 
n. However, the errors in the estimates increase with n as well. One must choose an optimal 
n such that f(z) is adequately parametrized, while the errors on the estimated amplitudes 
of f(z) are not too large to be useful. We've experimented with n > 10 parametrizations of 
f(z), and found that n = 10 is a good choice. 

The biased estimates of the dark energy density f q (z) and fk(z) in Fig.2 is mainly due 
to the bias in the estimate of Q m . Fig. 3 shows f q (z) and fk(z) estimated assuming that we 
know Q m = 0.3, with the same line types as in Fig.2. Even with this ideal assumption, it is 
only possible to marginally differentiate between the two models. 

Fig.4 shows the effect of adding a systematic shift of dm sys z to 100 random data sets 
with a realistic dispersion of 0.2 magnitudes. The line types are the same as in Fig.2. 
We have added a systematic shift in /x of (a) 0.01 z magnitudes; (b) 0.05 z magnitudes. 
Clearly, systematic shifts can significantly increase the bias in the estimate of Q m and the 
estimates of f(z) for z > 1.2, while having little effect on the estimates of f(z) at z < 1.2. 
Systematic errors as a function of z may arise from intrinsic properties of the supernovae 
varying with z such as progenitor evolution or dust characteristics changing with metalicity 
of the Universe. Accuracy may also be limited by the observations themselves in the case 
of k-corrections or selection biases. The estimated errors from these sources in the current 
surveys are between 2% and 10% ( [Schmidt et al. 1998 ). It is perhaps not realistic to expect 
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to measure both Q m and f(z) to high redshift accurately from SN la data alone. 

It is clear from Fig. 4 that even in the presence of plausible systematic effects, we 
can expect to measure the time- variation in the dark energy density f(z) with reasonable 
accuracy to a redshift of about 1.2. The bias and the errors in the estimated f(z) increase 
substantially beyond z = 1.2. 

We only applied our adaptive iteration method to 100 random samples, and with a 
resolution of AQ m = 0.02, because this method takes several hours per sample on a fast 
Sun work station in finding the best fit f(z) and Q m . It is work planned for the future to 
improve this promising method for application to much larger number of random samples, 
as well as adding fl x as an estimated parameter. 



6. Implication for dark energy models 

Recently, there has been a great deal of activity in exploring the possibilities of the 
existence of exotic dark energy ( Peebles fc Ratra 1988| , [Frieman et al. 1995| , |Caldwell, Dave 



fc Steinhardt 1998| , |Sahni fc Wang 2000| ) in the universe. While the present observational 
data are consistent with the dark energy being a cosmological constant, they do not rule 
out alternatives in the form of various scalar fields. 

It is important that we measure the time dependence of the dark energy density. If the 
dark energy density is measured to be constant in time within reasonable uncertainties, a 
cosmological constant should be favored, and more theoretical efforts should be directed 
toward the derivation of a cosmological constant from first principles. At the very least, 
this places strong constraint on the classes of scalar-field models for the dark energy. On 
the other hand, if the time-dependence of the dark energy density is established by the 
observational data, we would come to the exciting discovery of new physics in the universe. 

We have studied two dark energy models f g (z) and fk{z) (see Eq.(p2|)), representing 
two general classes of dark energy models, quintessence ( |Caldwell, Dave, fc Steinhardt 



19981) and k-essence (|Armcndariz-Picon, Mukhanov, fc Steinhardt 2000|) . This allows us 



to examine how well different models can be differentiated by realistic data, as well as the 
robust determination of the time- dependence of the dark energy density. 

Realistic future SN data, as described in the previous section, has the potential of 
determining the time-dependence of the dark energy density. This will clearly have a 
dramatic impact on models of the dark energy. 
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7. Conclusions 

To access the prospects of measuring the time variation in the equation of state, we 
have developed a promising adaptive iteration method that is powerful in extracting the 
dark energy density f(z) = px(z) / px(z — 0) from realistic data. This method is based on 
the requirement that the weak energy condition (energy density of matter is nonnegative 
for any observer) is satisfied. 

We have found that current type la supernova (SNe la) data are consistent with a 
cosmological constant, with large uncertainties at z > 0.5. We show that Q m (assuming 
a flat universe) and the dimensionless dark energy density f(z) = px(z)/px(z = 0) can 
be measured reasonably well to about z = 1.5 using type la supernova data from realistic 
future SN la pencil beam surveys, provided that the weak energy condition (energy density 
of matter is nonnegative for any observer) is imposed. For z > 1.5, the errors increase 
significantly, while the estimates become more biased, making it impossible to differentiate 
between different models (say, quintessence and k-essence). However, the correct trend in 
the time-dependence of the dark energy density can be clearly detected out to z — 2, even 
in the presence of plausible systematic effects. This would allow us to determine whether 
the dark energy is vacuum energy, or some exotic form of energy with a time-dependent 
density. 

The simulated data we used are for a SN pencil beam survey (either from ground or 
from space) with a square degree field of view and for an effective observational period of 
one year ([Wang 2000a| ), combined with SN data at smaller redshifts, so that there are a 
minimum of 50 SNe la per 0.1 redshift interval at any redshift. Although the dispersion of 
0.20 magnitudes (intrinsic plus observational) assumed in our simulated data is appropriate 
for ground-based surveys, we expect our results to apply qualitatively to space based 
SN pencil beam surveys as well, because our assumed dispersion of 0.20 magnitudes is 
dominated by intrinsic dispersion (about 0.17 magnitudes). 

At the completion of this lengthy numerical study of the feasibility of measuring 
the time-dependence of the dark energy density from realistic SN la data, we became 
aware of the recent paper by Maor, Brustein, & Steinhardt (2000). They claimed that 
distance-redshift relations derived from SNe la and similar classical measures are poor 
methods for resolving the time- dependence or measuring the amplitude of the equation 
of state wx(z) = px(z)/px( z ), and consequently no useful information can be obtained 
about the future fate of the universe. Our work has confirmed the difficulty of measuring 
the properties of the dark energy from realistic SN la data. However, we have found their 
assessment to be overly pessimistic. Their work indicates that it is impossible to tell 
whether the equation of state Wx(z) varies in time (see also, |Barger fc Marfatia 200 1|), but 



knowing that the dark energy density px{ z ) varies in time is sufficient to rule out vacuum 
energy as dark energy, thus giving support to exotic dark energy models. Our work has 
shown that one can indeed clearly detect the time-dependence of px (z) using realistic future 
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SN la data. 

The main problem reported by Moar et al. was the "smearing" effect of the 
multi-integral relation between the luminosity distance di{z) and the equation of state 
wx(z). Instead of di(z) and wx{z), our analysis uses the time derivative of the comoving 
distance r'(z) and the dimensionless dark energy density f(z). Thus our results are less 
affected by the smearing effect. Our method will be refined and made more efficient, and 
should become quite useful in analyzing future SN data. 

In view of our results, it is important that reasonable yet substantial observational 
efforts are devoted to future SN la surveys, for example, a dedicated 4m telescope for 
a SN pencil beam survey ( [Wang 2000a|) , combined with surveys of nearby SNe la. The 



total cost of such surveys would be modest compared to the great scientific return, the 
determination of the systematic uncertainties of SNe la as cosmological standard candles, 
and the measurement of the time- dependence in the dark energy density of the universe to 
constrain fundamental physics. 
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Fig. 1. — The dimensionless dark energy density f(z) (as parametrized by Eq.(|13D) measured 
from the current combined data of 92 SNe la. 

Fig. 2. — The measured dark energy density f(z) [see Eq.([T3|)1 in the absence of lensing 
noise and systematic shifts, when we apply our adaptive iteration method to 100 random 
data sets with a realistic dispersion of 0.2 magnitudes, for (a) n = 10, and (b) n = 6. The 
thick and thin solid lines are the assumed true f q (z) and fk(z) respectively. The circles and 
triangles are the estimated values of f q (z) and fk(z) respectively. 

Fig. 3. — The dark energy densities f q (z) and fk(z) estimated assuming that we know 
Q m = 0.3, with the same line types as in Fig.2. 

Fig. 4. — The effect of adding a systematic shift of dm sys z to 100 random data sets with 
a realistic dispersion of 0.2 magnitudes. The line types are the same as in Fig.2. We have 
added a systematic shift in fi of (a) 0.01 z magnitudes; (b) 0.05 z magnitudes. 
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Fig. 1. — The dimensionless dark energy density f(z) (as parametrized by Eq.(|l3D) measured 
from the current combined data of 92 SNe la. 
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Fig. 2. — The measured dark energy density f(z) [see Eq.flT3|)1 in the absence of lensing 
noise and systematic shifts, when we apply our adaptive iteration method to 100 random 
data sets with a realistic dispersion of 0.2 magnitudes. The thick and thin solid lines are 
the assumed true f q (z) and fk{z) respectively. The circles and triangles are the estimated 
values of f q (z) and fk(z) respectively, (a) n = 10. 
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Fig. 3. — The dark energy densities f q (z) and fk(z) estimated assuming that we know 
Q m = 0.3, with the same line types as in Fig.2. (a) n = 10. 
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Fig. 3. — (b) n = 6. 
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Fig. 4. — The effect of adding a systematic shift of dm sys z to 100 random data sets with a 
realistic dispersion of 0.2 magnitudes. The line types are the same as in Fig.2. (a) We have 
added a systematic shift in /i of 0.01 z magnitudes. 
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Fig. 4. — (b) We have added a systematic shift in /i of 0.05 z magnitudes. 



